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Abstract; Results of three-dimensional Monte Carlo calculation of low 
energy flux of atmospheric neutrinos are presented and compared with earlier 
one-dimensional calculations [1,2] valid at higher neutrino energies. These low 
energy neutrinos are the atmospheric background in searching for neutrinos from 
astrophysical sources. Primary cosmic rays produce the neutrino flux peaking at 
near E 40 MeV and neutrino intensity peaking near £7^=100 MeV. Because 
such neutrinos typically deviate by 20°~ 30° from the primary cosmic ray direc- 
tion, three-dimensional effects are important for the search of atmospheric neutri- 
nos. Nevertheless, the background of these atmospheric neutrinos is negligible for 
the detection of solar and supernova neutrinos. 

1. Introduction 

Recently one-demensional Monte Carlo calculations of cosmic ray produc- 
tion of neutrinos in the Earth’s atmosphere, including geomagnetic and solar 
modulation effects, were reported [1,2]. These calculations agree well with the 
flux and angular distribution of neutrinos of energy E„> 200 MeV observed in 
underground detectors. But these earlier one-dimensional cascade calculations are 
inapplicable to neutrinos of lower energy E „< 50 MeV important in neutrino 
astronomy. In fact, we shall see (Fig. 3) that the mean angular deviation 
between the primary cosmic rays and neutrinos is appreciable even for neutrinos 
of several hundred MeV. The flux of low-energy atmospheric neutrinos we obtain 
(Fig. 2) in the present three-dimensional calculation is negligible compared with 
the known flux of solar or supernova neutrinos, but maybe significant in the case 
of other sources of astrophysical neutrinos. 

The three-dimensional atmospheric neutrino flux is 

dN u {E v , 0 V> (f> v ) / dE „= J y„ U ( E p , 0 p , <f> p> X)(dN p /dE p )dE p du p) 

where y t ,(E J/ , 0 U , <t>„,E p , 0 p ) is the yield of neutrinos of energy E u , zenith angle 
9 V and azimuth angle <!> v by primary cosmic rays of energy E p and zenith angle 
<t > p . dN p /dE p is the primary cosmic ray spectrum, and fi ( E p ,0 p ,<f > p ) is the 
geomagnetic cut-off. This geomagnetic cut-off depends on geomagnetic latitude A 
and magnetic rigidity R = pc /e where p is the primary cosmic ray’s. , momentum. 
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2. Calculational details 

We modified the Gaisser-Protheroe-Stanev one-dimensional hadron interac- 
tion model [5] by assigning to secondary particles the transverse momentum dis- 
tribution. 

W ( a,x t )=(a +l)(a +2)x t (l+x t ) a . 

Here a =2 p /< p p > -3 for p > 1.5< p t > , a =0 for p < 1.5< p t > , x t -p t /p and p 
is the incident hadron’s momentum. At high p, W ~ p t exp(—ap t ). The angular 
deviation of secondary hadrons from the incident hadron direction produced by 
the above formula affects our results insignificantly. 

In our computer program analytic formulae are used for the energy distribu- 
tion of decay particles. For those energies, a microcanonical ensemble average is 
taken; total energy is conserved in each individual decay. The secondary direc- 
tions are, however, assigned according to a canonical ensemble; the total momen- 
tum is conserved on average but not in individual decays. Interaction cases are 
treated similarly. We use energy-independent energy loss rates for charged parti- 
cles by air ionization in a simple isothermal atmosphere to determine the decay 
height and energy, and energy-dependent formulae to determine the interaction 
height and energy. Our results are insensitive to these parameters. 

For the comparison with the earlier one-dimensional calculation we calcu- 
lated the yield function 

¥* 

Y i/{E & p ) = J' V i/{E 0 v> fii/jEp } Op )d Up,, 

integrated over appropriate neutrino directions. 

3. Results and Conclusions 

Fig.l shows the yield Y v of vertically incident protons of several energies 
E p , integrated over the downward 2n solid angle of neutrino’s directions. 
Because low-energy neutrinos are made mainly from low-enrgy pions and muons 
at rest, for all primary cosmic ray energies the peak flux occurs near 40 MeV and 
the peak intensity near 100 MeV. Because most muons decay at rest and 
< E v > =0.3 E p,, < Ep, > =0.35 E M , the atmospheric neutrino flux peaks near 40 
MeV.' 

Fig.2‘shows the downward neutrino flux produced by vertically incident pri- 
maries at high geomagnetic latitude, where 

dNjdE* J Y u n(dN p /dE p )dE p . 

While this expression is not exact for a three-dimensional cascade, it compares 
simply with the one-dimensional calculations. The near-vertical neutrino flux is 
also insensitive to zenith angle. Because, if not supressed by geomagnetic cut-off, 
the primary cosmic rays have a steep power-law spectrum, the main contribution 
to the low-energy neutrino flux is from primary cosmic rays just above the pion- 
production threshold energy. The total flux is therefore very sensitive to geomag- 
netic cut-off [6]. 



Flux dN v jdE u (cm 2 sr s GeV) 1 Yield /(p ro ton.GeV) 
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Fig. 1. Integrated neutrino yields from vertically 
incident protons. 2 5 10 20 50 GeV cases from 

bottom to top. 
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Fig. 2. Calculated downward neutrino flux at high 
geomagnetic latitude 
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For neutrinos above E^—200 MeV, this three-dimensional calculation agrees 
with earlier one-dimensional calculations [1,2]. At lower energies ( E v < 50 MeV), 
however, we now obtain neutrino fluxes significantly lower than those one- 
dimensional calculations would give. 

This atmospheric neutrino background is totally negligible compared to the 
solar neutrino flux and the flux of neutrinos expected from supernova at any rea- 
sonable distance (< 10 3 Mpc). 

Fig. 3 shows the average angular deviation of down-going u e s of energy E v 
produced by vertically incident protons of energy E p =2, 10, 50 GeV. Other neu- 
trino types have similar spectra. While average angular deviation is 20° ~ 30°, 
neutrino direction at neutrino energy below 50 MeV is near isotropic and a long 
tail of more energetic neutrinos deviating by 5°~ 10° is produced at E v > 200 



Fig. 3. i/ e angular deviation from primary proton. 

2 GeV (O), 10 GeV (A), 50 GeV (+) cases. 
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